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Using the "Quality Factor" (QF) method, we analyse the scaling properties of deep-inelastic pro- 
cesses at HERA and fixed target experiments for x < 1CP 2 . We look for scaling formulae of the 
form <r 7 p (t), where t(L = log Q 2 , Y) is a scaling variable suggested by the asymptotic properties 
of QCD evolution equations with rapidity Y. We consider four cases: "Fixed Coupling", corre- 
sponding to the original geometric scaling proposal and motivated by the asymptotic properties of 
the Balitsky-Kovchegov (BK) equation with fixed QCD coupling constant, two versions "Running 
Coupling 1,11" of the scaling suggested by the BK equation with running coupling, and "Diffusive 
Scaling" suggested by the QCD evolution equation with Pomeron loops. The Quality Factors, quan- 
tifying the phenomenological validity of the candidate scaling variables, are fitted on the total and 
CNj , DVCS cross-section data from HERA and predictions are made for the elastic vector-meson and for 

. the diffractive cross-sections at fixed small xp or (3. The first three scaling formulae have compa- 

rably good QF while the fourth one is disfavored. Adjusting initial conditions gives a significant 

| improvement of the "Running Coupling II" scaling. 

*so : 

I. INTRODUCTION 

Geometric scaling l] is a remarkable empirical property verified by data on high energy deep inelastic scattering 
(DIS) i.e. virtual photon-proton cross-sections. It was realized that one can represent with reasonable accuracy the 
'— cross section cr 7 p by the formula a 1 P (Y, Q) — a 7 (r) , where Q is the virtuality of the photon, Y the total rapidity 
' in the 7*-proton system and 

T = logQ 2 -logQ2(y)=logQ 2 -Ay , (1) 

00 ' 



is the scaling variable. The value A ~ 0.3 confirms the value found within the Golec-Biernat and Wusthoff model 
£Nj ' where geometric scaling was explicitely used for the parametrization. 

| The scaling using the variable r defined in Formula (|TJ) was intimately related to the concept of saturation 0], i.e. 
the behaviour of perturbative QCD amplitudes when the density of partons becomes high enough to exercise nonlinear 
effects ensuring the unitarity bound. Indeed, there were many theoretical arguments to infer that in a domain in Y 
and Q 2 where saturation effects set in, geometric scaling is expected to occur. Within this framework, the function 
Q S (Y) can be called the saturation scale, since it delineates the approximate lower bound of the saturation domain. 
• i-h , This type of geometric scaling is motivated by asymptotic properties of QCD evolution equations with rapidity. 
Using the nonlinear Balitsky-Kovchegov (BK) equation [3| which represents the "mean-field" approximation of high 
energy (or high density) QCD, geometric scaling could be derived from its asymptotic solutions Q. This equation is 
supposed to capture some essential features of saturation effects. Considering the BK equation with fixed coupling 
constant leads asymptotically to the original geometric scaling of Formula {1} . Considering a running coupling leads 
to the following scaling: 

r = log Q 2 - log Q 2 S (Y) = log Q 2 - \W , (2) 
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Recently [6|, it was noticed that the scaling solution ^ of the BK equation with running coupling is only approx- 
imate and not unique. Another equivalent approximation leads to a different scaling variable, namely 



r = logQ 2 -A — ^— . (3) 
logQ 2 

The key ingredient to theoretically prove geometric scaling of the asymptotic solutions of the nonlinear BK equation 
is the traveling wave method [5(. Indeed, the BK equation admits solutions in the form of traveling waves M{L—v c t). 
L = log Q 2 has the interpretation of a space variable while t, interpreted as time, is an increasing function of rapidity 
Y, namely tec Y for the fixed coupling case and t oc VY for the running one. v c is the critical velocity of the wave. 
These properties confirmed results previously obtained Q, for instance by replacing the nonlinear damping term by 
absorbing boundary conditions, and thus considering only the linear part of the BK equation, which is equivalent to 
the Balitsky Fadin Kuraev Lipatov (BFKL) equation [8(. 

The subsequent elaboration on QCD evolution equations led to go beyond the mean-field approximation. The effect 
of fluctuations was examined in Ref. 0] in the fixed coupling scheme and gives rise to a new "diffusive scaling" , the 
scaling variable being 

-- l0gQ2 - Ar . (4) 



Y 

The aim of this paper from a theoretical point-of-view is to test and compare the different scaling behaviors, arising 
from different versions of QCD evolution, using the data available from HERA. We shall study the phenomenological 
relevance of the four kinds of scaling and refer to them in the following as "Fixed Coupling" for the variable |T]), 
"Running Coupling I" and "Running Coupling II" for the variables ([2]) and |3]) respectively, and to "Diffusive Scaling" 
for Our method can be easily applied to any new proposal of scaling. 

From an experimental point-of view, one wants to include not only the newest published data on a 1 p from HERA, 
but also from other processes where the small- a; physics allow one to discuss scaling properties. We thus include 
the data on DVCS cross sections and, following the geometric scaling analysis performed in Ref. we extend the 
analysis to elastic vector- meson production and to diffractive cross-sections at fixed and small xp. Moreover, and 
as a new aspect of scaling analysis, we are led to also include diffractive cross-sections at fixed and small (3, since 
scaling should also be considered for the 7* Pomeron cross section, following the Ingelman-Schlein interpretation of 
diffractive structure functions [ill ]. 

In order to define how good the scalings are, the authors of Ref. [ll| introduce a Quality Factor (QF) as an 
estimator on the validity of scaling. The main property of the QFs is that it does not depend a priori on a given 
parametrization of the scaling curve as in an ordinary x 2 approach. The aim of this paper is thus to extend on 
the same footing the QF studies to all scaling laws, including the new scaling Running Coupling II, and to compare 
them. We perform a systematic analysis using the QF techique investigating the scaling properties in deep-inelastic 
scattering using all data available. It will include both the largest available set of published data where scaling is 
expected to occur and the new status on theoretical scaling properties using QCD evolution equations. 

The plan of the paper is as follows. In section II we describe the theoretical motivation and the give the precise 
formulation of the various scaling hypotheses. In section III we present the Quality Factor method we shall use in 
the analysis and the set of data we consider for the fitting procedure. We provide and comment our results in section 
IV, and present our conclusions in the final section V. 



II. SCALING VARIABLES 



Let us sketch the theoretical motivation for the different forms of scaling HQ in deep-inelastic scattering and their 
extension to DVCS cross-section data, elastic vector-meson and for diffractive cross-sections at fixed and small xp or 
(5. We shall focus on the general arguments leading to scaling which could be independent from too specific theoretical 
predictions. The different types of scaling we refer to originate from properties of nonlinear QCD evolution equations, 
and the existence of asymptotic traveling wave solutions. We shall first introduce the general hint for these traveling 
wave solutions in the four cases and subsequently study the effect of varying the initial values of the phase space 
variables Y and L. This amounts to introduce some natural non asymptotic effects. 

For the three first cases (0 [3 HI , let us start with a "mean-field" equation for the dipole-target amplitude of the 
type 

d Y N{L,Y) =a X {-d L )N{L,Y)-aN 2 {L,Y), (5) 

where L — \og{k 2 / k 2 ) , k is the gluon transverse momentum (Fourier conjugate to the dipole size), kg is an arbitrary 
constant and Y is the rapidity. In Eq. |1J , the coupling constant will be considered in the following as fixed, or running 
such that a(L) — 1/bL , b = llN c — 2Nf /12N C where Nf and Nc are the numbers of flavours and colours respectively. 
As we will show, for the sake of a general derivation of scaling, the differential operator kernel x(— can be considered 
to be general, provided the function x(l)/l admits a minimum value at some point j c . Indeed, the kernel at the leading 
logarithm (LL) accuracy of the perturbative QCD expansion, the BFKL kernel [8] x(7) = 2-0(1)— -0(7)— -0(1— 7) verifies 
this propert y b ut large corrections at next-to-leading logarithm (NLL) accuracy exist (see, for the initial references, 
e.g. [HI, dHHH, [IB]). We will thus consider a general argument, instead of using a specific form of the QCD kernel. 

Let us first consider the "Fixed Coupling" scheme a = est. The solution to the linear part of equation §E§ corresponds 
to a linear superposition of waves: 



^F) = / c |l^(7)e X p{-7(i:-^y 



(6) 



where 0^(7) = ctxi'j) is the Mellin transform of the kernel. In particular, one can interpret §E§ by stipulating that each 
partial wave of wave-number 7 has a velocity ^(7) = By contrast, the critical velocity is defined @, [l3| by the 



minimum j c of the exponential phase factor v 



dry 



;(7e) 

7c 



v c . v c is thus the minimal velocity of the waves. 



So, the balance between the initial velocity (for fast enough initial condition with 70 > 7 C , see [f|) and the damping 
due to the nonlinear term in Eq ^ leads to an asymptotic blocking of the velocity (and of the wave front) on the 
minimal and critical value v c . As a consequence, reporting this dominant value in formula (j6]), one finds 

W,Y)~e X J- 7c (L-^Y)} , (7) 



which is the expression of scaling with the variable |T]). 

For the running coupling case of Eq. 5, the situation is more subtle, since different approximate traveling wave 
solutions may coexist. In the original derivation of Ref. one starts with the linear solution analogue to Eq. ((U, 
namely 



where X(j)(= J 1 drf xil')) pl avs the role of the new effective kernel. The key point is that the equivalent of the 
time variable is now y/Y instead of Y. This VY comes from a saddle-point integration over w, the conjugate variable 
to the rapidity (l5j . Consequently, after the same "blocking" mechanism due to the nonlinear damping terms, the 
scaling variable is t(L, Y) = logQ 2 — X\/Y , which is nothing else than the "Running Coupling I" case, see Eq. 
However, the scaling inferred by Eq. {5]) is neither an exact nor unique approximate scaling solution. 

As an alternative to that Mellin transform method, it was proposed recently in Ref. 0] to search directly for the 
running coupling analog of partial wave solutions, i.e. scaling solutions of Eq. ([5]) with running coupling a(L) — 1/bL, 
with some generic velocity v. Imposing the scaling of the left hand and right hand side of the equation as a function of 
the scaling variable t(L, Y) leads to two different constraints, namely the functions 6L9yA/"[r(L, Y)] and <9£A/"[t(L, Y)] 
should be both functions of r only 

bLd Y U[r(L,Y)} = f(r) (9) 
dlN[r{L.Y)] = g(r) (10) 



which are conditions to have a scaling solution Af(L,Y) = J\f[r(L,Y)]. In fact conditions © and (|10p cannot be 
fulfiled simultaneously (contrary to the 'Fixed Coupling" case, see (6() and can only be considered in an approximate 
way. If one choses to fulfil condition (fit)]) exactly, and © only approximately, one finds again the scaling variable 
t(L,Y) = logQ 2 — XVY , which is the "Running Coupling I" case, see Eq. ((2)). By this way one recovers the scaling 
variable previously found in Refs. [H, 0] in a different way. If, conversely, one chooses to satisfy exactly (jH]) and only 
approximately (|10p . one finds a new form of scaling, with the variable t(L,Y) = L — , which corresponds to 
the "Running Coupling II" case, see Eq.©. As the "blocking" mechanism is more general than the approximations 
leading to the variables "Running Coupling I" or "Running Coupling II" , one expects that Eq. (|5|) with running 
coupling has no exact scaling solution, but admits at least two different approximate solutions, namely "Running 
coupling I" and "Running coupling II" . One aim of our study is to make a phenomenological comparison between 
both solutions. 

For the case of diffusive scaling, one considers the modification of Eq. ([5]) due to Pomeron-loop contributions, which 
appear when the system of QCD partons or dipoles is dilute. It may seem irrelevant to an equation describing the 
high density partonic effects on the amplitude, but the structure of the equation, with an exponential increase of the 
linear regime may lead to important modifications 1 w.r.t. the mean-field solutions of Indeed, it was shown that, 
while traveling-wave solutions of the type (O are formed, the effect of event-by-event fluctuations leads to a stochastic 
superposition of these waves around an average solution which is governed by some diffusion coefficient and thus to 

a different scaling variable t(L,Y) = log;( ^-I AY which corresponds to the "Diffusive Scaling" case, see Eq. ([!]). 

The following comments are in order. One would like to confront the QCD predictions on the scaling form of the 
amplitude (i.e. the wave front) with data. However, the precise determination of the wave front seems at present 
to be less reliable than the scaling variable (i.e. the wave structure). The form of the linear kernel and of the non 
linear damping at next-to- leading logarithmic order have just been derived . However, in the kinematical range 
where scaling is observed, higher order contributions may be relevant, leading to a modified kernel in Eq. (5). Hence 
we will stick to a rather general and model independent (for each case) prediction for the scaling variable. The 
phenomenological analysis may on contrary help the development of the theoretical investigation. 

Another important comment concerns the non-asymptotic corrections. While some of them can be deduced from 
the blocking mechanism itself (see e.g. [j|) we shall consider the simplest and expected nonasymptotic corrections 
due to the definition of the typical momentum and rapidity scales A and Y . Hence we will consider, if necessary 
the modification of all scaling variables JIE]), by a shift L — > L — log A 2 and Y — Yq. Note that the reference value 
(without shift) is always Yo = 0. 

Let us now consider diffractive, elastic vector-meson and DVCS cross-sections. As discussed in Ref. [lOj, geometric 
scaling was shown to be valid and the "Fixed Coupling" was verified within uncertainties for these reactions, using 
xp instead of the Bjorken x in the scaling variable fl}. 

In our phenomenological discussion, we shall thus extend the scaling properties to all four cases including the one 
studied in [l(| , by considering the various scaling variables (PQl . 

We shall thus investigate the following scaling properties. The scaling prediction is 

°iZc7(^ Q 2 )= a Dvcs p ( T i x > Q 2 \) (ii) 

for the DVCS cross-section, 

4ff V *l?r> Q 2 ' M v) = Vv£T Vp (tv[xp,Q 2 + My])) (12) 

for the elastic vector meson cross sections, where the hard scale is assumed to be Q 2 + M V , from the known properties 
of the vector meson wave function, and 

— ^ ((3,x F ,Q 2 ) = ^l (P,T d [x F ,Q 2 }) (13) 

for the diffractive cross-section at fixed j3 and small xp. 

For completion, besides the cases studied in [13], we may also consider the diffractive cross section at fixed xp and 
small [3, namely 

-^LL (f3, X p,Q 2 ) = -ML (xp,r b [/3,Q 2 }) , (14) 



1 Recent results on a toy model with running coupling and fluctuations [lgj seem not to give diffusive scaling. A more complete QCD 
study is still missing. 





scaling 


t formula 


parameters 


FC 


"Fixed Coupling" 


log Q 2 - XY 


A 


RCI 


"Running Coupling I" 


\ogQ 2 -XVY 


A 


RCIb 


"Running Coupling Ibis" 


log<y - XVY-Yq 


A, Y 


RCII 


"Running Coupling II" 


log(Q 2 /A 2 ) - Aj^^y 
A =0.2 GeV 


X 


RCIIb 


"Running Coupling Ilbis" 




A, Y , A 


DS 


"Diffusive Scaling" 


log(Q z /A^)-\Y 

A =1 GeV 


A 


DSb 


"Diffusive Scaling bis" 




A, Y a , A 



TABLE I: Scaling variables used in the fits to deep inelastic scattering data. 



since in the hypothesis of factorisation of the Pomeron flux, it relates to c 7 v ^ xv (Ji, Q 2 ) (i.e. the total "Pomeron- 
photon" cross-section) at small (3 for which the scaling arguments should also hold. 

The prediction for DVCS (see Eq. (jlip) does not rely on any nonperturbative assumption. This interesting feature of 
DVCS cross-sections will allow us to include the data in the fitting procedure. Note that we will not use the diffractive 
data nor the vector meson ones for the fit. Instead, we will consider the fit to total and DVCS cross-sections as scaling 
predictions for diffractive and vector meson data. 

Hence, in our phenomenological discussion, we shall extend the scaling properties to our four cases, by considering 
the various scaling variables (OH in the equations (|1H114|) . 



III. THE QUALITY FACTOR METHOD 



In order to fit the parameters in the scaling variables ([IH]), it is useful to use a quantity called quality factor (121 ] 
(QF). The quality factor is used to find the best parameters for a given scaling and to compare quantitatively the 
different scalings. 

Given a set of data points (Q 2 ,x,a = a(Q 2 ,x)) and a parametric scaling variable r = t(Q 2 ,Y — log(l/x);X) 
we want to know whether the cross-section can be parametrised as a function of the variable r only. Since we do 
not know the function of the r variable describing the data, we need to define the QF independently of the form of 
that function. Therefore the QF is defined in such a way that it quantitatively describes how close the data points 
a = a(Q 2 ,x) are to the scaling law er = <t(t), for a given r = t(Q 2 ,Y; A). 

Let us consider a set of points (ui, Vi), where Mi's are ordered, and introduce the QF as follows [HJ 



QF(X) = 



E 



(Ui - Ui-i) 2 + < 



(15) 



where e is a small constant that prevents the sum from being infinite in case of two points having the same value of 
u. Using this definition, the contribution to the sum in (|15j) is large when two successive points are close in u and far 
in v. Thus, we expect a set of points lying close to a unique curve to have large QF (small sum in (TL3)) ) compared to 
a QF of a set of points that are more scattered. 

Since the cross-section in data differs by orders of magnitude and r is more or less linear in log(Q 2 ) (see HE]), we 
decided to take Vi = log(ai) and Ui = t^(A). This ensures that the low Q 2 data points contribute to the QF with a 
similar weight as the higher Q 2 data points. To complete the definition of the QF, the set (ui,V{) is also rescaled so 
that < Ui,Vi < 1, and ordered in u before entering the QF formula. All the QFs in this paper are calculated with 
e = 0.01. 

In order to test a scaling law r, we search for the parameter A that minimises the 1/QF variable. To do that, we 
fit 1/QF with the MINUIT package. Given the maximum value of the QF, we are able to directly compare different 
scaling laws. In this paper, we test the scaling laws (OH). The different scaling variables are given in Table |U 

Some comments are in order. We always define two versions of the scaling variables, the former depending only on 
a single parameter A, and the later on two or three parameters A, Yq, A. As an example, we define in this way the 
"Running Coupling II" and "Running Coupling Ilbis" in this way respectively. The "Fixed Coupling" scaling does 



not depend on the scale A and the rapidity Yq since these additional variables disappear in the definition of the quality 
factor QF. The conclusion remains the same for the A dependence of the "Running Coupling I" scaling. Concerning 
the "Running Coupling II" scaling, we fixed the parameter A to a typical non perturbative scale inside the proton, 
A = A Q cd = 0.2 GeV. 

The QF defined above raise two additional technical issues which require further checks. First, it does not take into 
account the errors on data points. The errors can be introduced to the QF in an error weight function by modifying 
the definition (|15|) in the following way 



QF{\) = 



E 



Ui-i) 2 W M _i 



i -i 



(Ui - Ui-l) 2 + £ 2 



(16) 



where Wjy is the error weight function of the data points In the simplest form, we can take 



Wij = fnpj, (17) 

where pi is the relative statistical or uncorrelated error of point i. This definition allows to take into account the 
statistical scatter of the different data points. This is particularly important if the statistical uncertainty differs from 
one data point to the other, and is significantly large. The data points we use in the following (mainly the structure 
function F2 from the HI and ZEUS experiments) show very little statistical uncertainties. We compared the results 
using either formula (| 1 5 j) or formula (|16|). and as expected, we get similar values of the parameters for both definitions. 
For sake of simplicity, we choose the QF ([15)) without any error weight function. 

The second possible shortcoming is an ambiguity of the QF related to the (ut, Vi) point ordering if two or more Uj's 
are equal. The order of these points in the data set is not unique; any permutation of these points satisfy the definition 
above. Obviously, we get diffrent QF for different permutations. The comparison of the fit results with different data 
point ordering was performed and the results were found similar (both in the value of the fitted parameters and the 
QF itself). 



IV. FITS TO F 2 AND DVCS DATA 

The first natural data set to test the different scaling laws are the F2 data measured at low x by the HI [13] , ZEUS 
[2lj | . NMC [13] and E665 [23| experiments. These data are very precise and cover a wide range in x and Q 2 . 

We choose to consider all available data from HI, ZEUS, NMC and E665 experiments with Q 2 in the range 
[1;150] GeV 2 and x < 0.01. Fig. [T] displays the different kinematical domains reached by the HI, ZEUS, NMC and 
E665 experiments. While the HI and ZEUS kinematic coverage is similar (with a tendency for HI to reach lower x 
and Q 2 values), the data from NMC and E665 cover mostly the low Q 2 , higher x kinematic domain of our sample. We 
choose not to consider the data with x > 10~ 2 since they are dominated by the valence and sea quark densities, and 
the formalism of saturation does not apply in this kinematical region. In the same way, the upper Q 2 cut is introduced 
while the lower Q 2 cut ensures that we stay away from the soft QCD domain. We will show in the following that 
the data points with Q 2 < 1 GeV 2 spoil the fit stability. Two kinds of fits to the scaling laws are performed, either 
in the full mentioned Q 2 range, or in a tighter Q 2 range [3; 150] GeV 2 to ensure that we are in the domain where 
perturbative QCD applies. We will study the Q 2 -dependence of the fit parameters, as a possible indication of smooth 
scaling violations. 



A. Fits to F 2 data with Q 2 > 3 GeV 2 



We first give the results of the different scalings for the structure function F2 with Q 2 > 3 GeV 2 . Fig. [2][6] show the 
normalised QF dependence on the value of the fitted parameters and the scaling curves in the Q 2 > 3 GeV 2 range 
for the "Fixed Coupling" , "Running Coupling I" , "Running Coupling II" , "Running Coupling II bis" , and "Diffusive 
Scaling" respectively. To see the impact of the different data sets, we choose to divide the data points in HI alone, 
ZEUS alone, Hl+ZEUS, and all data. In the figures, we normalise the QF to 1. so that we are able to compare 
the values of the A, Yq and A parameters for the different data sets. Table [H] give the fit parameters and absolute 
QF values — to distinguish the quality of the different scalings — for the various scaling laws (see Table I). On the 
different figures, we see that the A parameter is well determined in the fit and the different data sets lead to values of 
A which are quite close to one another. As an example, for fixed coupling, A ranges between 0.33 to 0.38, for "Running 




FIG. 1: Kinematical domain in x and Q 2 covered by the different data sets from HI, ZEUS, NMC and E665 used in this anlysis. 
Only the kinematical domain with x < 10~ 2 and 1 < Q 2 < 150 Gev 2 is shown. 



Coupling I" between 1.74 and 1.84, for "Running Coupling II" between 3.23 and 3.44, and for "Diffusive Scaling" 
between 0.31 and 0.37. The scaling plots (always on the right of FiglH IH andEJ show similar behaviour for the 
different scalings while it is clear that the scaling plot for diffusive scaling is slightly worse. 

To improve the scaling quality, we checked the dependence of our results to additional parameters which can be 
introduced in the scaling. Introducing a shift in rapidity Fj or a scale A in the "Fixed Coupling" scaling does 
not change the results since these new parameters disappear in the definition of QF. This effect is the same if one 
introduces a A scale in the "Running Coupling I" scaling. Introducing a new parameter lo to shift the rapidity in 
the "Running Coupling I" scaling does not lead to an improvement of the scaling. The fit leads to a value of Y 
compatible with for all data sets considered. The results are similar when one introduces the Yq parameter and 
the scale A in "Diffusive Scaling" and the parameters are found to be close to and 1 after the fit. We notice that 
"Diffusive Scaling" including the additional parameters Yq and A admits also other solutions leading to better QF, 
but these solutions are either instable from one data set to another or lead to unphysical values of the parameters 
(large values of A and large negative values of Yb). 

The QF plots in figure [5] show the dependence on one of the parameters while the remaining parameters are fixed 




FIG. 2: F2 data: Normalised QF to 1. as a function of A and scaling curve with A fixed to the value corresponding to the best 
QF for "Fixed Coupling" . A Q 2 > 3 cut was applied to the data. 



to the fitted values. It is important to notice that in the case of the "Running Coupling Ilbis" version, it is difficult 
to find the absolute minimum since the QF = QF(X, Yq, A) is not smooth and shows several minima. There is even 
a plateau of minimal values in A, Yq parameters and we picked the parameter values leading to the smallest values of 
QF. This is why the "Running Coupling Ilbis" scaling law shows larger spread concerning the parameter values (A 
varies between 3.91 and 4.89). The typical values for the lo (respectively A) parameters vary between -1.2 and -2.5 
(respectively 0.2 and 0.35). The values of A are in a typical range of non-pertubative scale in the proton as expected. 

We also notice in Table [ilj that the best quality factor is obtained for "Running Coupling Ilbis" for all data sets 
while the worst one is for the diffusive scaling. The comparison of the different QF is also shown in Fig. left. 



B. Fits to F 2 data with Q 2 > 1 GeV 2 

Figure [7] shows the QF plots for different scalings in the Q 2 range [1; 150] GeV 2 . We notice that the spread on the 
parameter values is larger than after the cut on Q 2 > 3 GeV 2 . The values of the parameters are given in Table [Til 
We also notice a tendency to find lower values of A for the "Fixed Coupling" , "Running Coupling I" and "Running 
Coupling II" scalings when one introduced more low Q 2 data. The QF values are similar for the "Fixed Coupling" , 
"Running Coupling I" , and "Running Coupling Ilbis" — with a tendency to be slightly better for "Running Coupling 
Ilbis" — and is worse for diffusive scaling as before. The corresponding scaling curves are displayed in Fig. [HI They 
show similar behaviour with a tendency for "Diffusive Scaling" to show a larger spread of data. These figures also 
depict the Q 2 < 1 GeV 2 data points in grey, which were not used in the fit. These data points tend to show the same 
kind of scalings with a larger spread, and this is why they were not included in the fit. 

Fig. [9J right, gives the QF = QF(X) curves for all data and different scaling laws in the Q 2 regions [1; 150] GeV 2 
(on Fig. [5] left in the [3; 150] GeV 2 range) so that the different scalings can be compared easily. As we mentionned 
already, in both Q 2 ranges, the "Fixed Coupling" , the "Running Coupling I" , and the "Running Coupling II" are 
quite good, and the "Diffusive Scaling" is disfavoured. For Q 2 > 3 GeV 2 , the "Running Coupling Ilbis" scaling shows 
the best QF. 




FIG. 3: F2 data: Normalised QF to 1. as a function of A and scaling curve with A fixed to the value corresponding to the best 
QF for "Running Coupling I" . A Q 2 > 3 cut was applied to the data. 



C. Dependence of the fitted parameters on Q 2 

To study the Q 2 dependence of the fitted parameters and to quantify the amplitude of scaling violations, the data 
points are divided into four separate Q 2 samples: [1;3], [3; 10], [10; 35], and [35; 150] GeV 2 . The number of points 
in each region is similar (see table HTH . The lower and upper bounds of 1 GeV 2 and 150 GeV 2 are chosen for the 
same reason as before. The Q 2 dependence of the parameter A is depicted in figure I1Q[ and the fitted parameters 
together with the QF values are shown in table IIIII There is a slight increase of the A parameter in the case of "Fixed 
Coupling" while "Running Coupling I" is quite flat. This can be easily understood since "Running Coupling I" shows 
a natural Q 2 evolution. We notice a stronger increase in the case of "Running Coupling II" , showing the breaking of 
this scaling as a function of Q 2 . Since this scaling gives already the best QF, it would be worth to study the breaking 
of scaling and introduce it in the model to improve further the description of the data. The A parameter decreases 
strongly (especially in the last Q 2 bin) for the diffusive scaling. The fact that A depends strongly on Q 2 confirms the 
fact that this scaling leads to the worst description of the data. 

D. Fits to DVCS data from HI and ZEUS 

After fitting all HI and ZEUS F2 data, it is worth studying whether the DVCS data measured by the same 
experiments [24( lead to the same results. The amount of data is smaller (34 points for HI and ZEUS requiring 
x < 0.01 as for F2 data) and the precision on the A parameter will be weaker. The kinematic coverage of the DVCS 
data covers a smaller region in x and Q 2 as F2: 4 < Q 2 < 25 GeV 2 and 5 • 10~ 4 < x < 5 • 10~ 3 . The results of the 
fits can be found in table IIVI and Fig[Tl] and [TSJ The scaling results are displayed only for "Fixed Coupling" and 
"Running Coupling Ilbis" since all plots look similar. To facilitate the comparison between the results of the fits to F2 
and DVCS data, a star is put in Fig. [12] at the position of the A value fitted to the Hl+ZEUS F2 data with Q 2 in the 
range [3; 150] GeV 2 . We note that the DVCS data lead to similar A values to the F2 data, showing the consistency of 
the scalings. The values of the QF show a tendency to favour fixed coupling, but all different scalings (even "Diffusive 
Scaling") lead to reasonable values of QF. 




FIG. 4: F2 data: Normalised QF to 1. as a function of A and scaling curve with A fixed to the value corresponding to the best 
QF for "Running Coupling II" . A Q 2 > 3 cut was applied to the data. 



V. IMPLICATIONS FOR DIFFRACTION AND VECTOR MESONS 



In this section, we check if the scalings found in the previous section can also describe diffractive and vector meson 
data. Since these data are much less precise than the Fi or DVCS data and depend more on non-perturbative inputs 
(meson wave function, diffractive parton distribution inputs...), we choose to impose the same values of parameters 
found in the previous section and check if the scaling is also observed using this value. A ht to the diffractive data leads 
to large uncertainties on the parameter values since the number of data points is quite small and their uncertainties 
large. 

Let us first describe the diffractive data. We use the latest t-integrated diffractive cross section measurements from 
HI and ZEUS experiment [25|, derived from the diffractive processes e+p^e + X + Y, where the proton stays 
either intact or turns into a low mass state Y. The diffractive cross section in terms of the ^-integrated reduced cross 
section a^^ 3 \xp,x,Q 2 ) reads 



d 3 a ep^eXY 



l-y+ V ^)<r? {3) {x P ,x,tf), (IS) 



dxpdxdQ 2 xQ 4 \ 2 

a D(3) _ pD(S) t_ pD(3) (1Q) 

" *2 l + (l_y)2^ • W 

The latter equation can be simplified to — F^ 3 \ since it is a very good approximation anywhere but at 

large y. HI large rapidity gap measurements are realised with My < 1-6 GeV cut, whereas ZEUS measures with 
the My < 2.3 GeV cut using the M^-method. Both experiments do not measure exactly the same cross section. 
However, the difference is a known constant factor; the ZEUS data points can be converted to the same My range 
as HI by multiplying ZEUS values by factor of 0.85. In case of the e + p — ► e + X + p processes where the proton is 
tagged in Roman pot detectors (both HI and ZEUS), the data points have to be multiplied by 1.23. 

In order to test the geometrical scaling properties we plot the r dependence of (3da^ Li JJ* Xv / 'df3, which can be 



expressed in terms of diffractive structure function as follows 



^ d/3 " "Q^^ 2 ■ (20) 



The definition of the r variable is the same as in Table I replacing the variable x by xp, Q 2 remaining the same. 

In figure [TBI we show the HI and ZEUS diffractive data for ftda'J li jJ* Xp / dft as a function of the fixed coupling 
scaling variable r for six different fixed values of ft: 0.04, 0.1, 0.2, 0.4, 0.65, and 0.90. Only the data points in the Q 2 
region [5; 90] GeV 2 and with xp < 0.01 arc plotted. The restriction to the mentioned fixed values implies that the 
data set is dominated by HI data. We just give in Fig. [T31 the results for "Fixed Coupling" since the other scalings 
give similar results which are not distinguishable given the large uncertainties on F 2 D . The parameters fitted to the 
HI and ZEUS F^ data are used to make the prediction on . Up to ft = 0.65 (the first five bins), the data points 
lie on a unique curve and thus confirm the geometrical scaling prediction. The data points at the highest ft values 
(low masses) are more sensitive to low mass resonances and are not expected to show a perfect scaling. The values of 
the QF are also given in Table fVl In many bins, "Running Coupling Ilbis" gives the best QF. 

Following the Ingclman-Schlcin model of the Pomeron, we also tested the scaling in ft, Q 2 in the diffraction data. 
The definition of the r variable is the same as in Table I replacing the variable x by ft, Q 2 remaining the same. Figure 
[Ml presents the HI and ZEUS diffractive data in the same Q 2 region with ft < 0.5 at five different fixed values of 
xp: 0.0003, 0.001, 0.003, 0.01, and 0.03. The restriction to these fixed values results in favouring the HI LRG data. 
Similarly as in the previous case, the parameter values from the fit to HI and ZEUS F2 data are used. The scaling is 
definitely not as good as the one for a fixed (3. 

The fits to HI and ZEUS data are also tested on vector meson data from both experiments (26[. Note the scale 
My present in the scaling function in (fl"2")) . The definition of the r variable is the same as in Table I replacing the 
variable Q 2 by Q 2 + M v My being the mass of the vector meson, x remaining the same. The QFs are given in table 
IVII and scaling curves for "Fixed Coupling" are given in Fig. [I5l We only show the results for "Fixed Coupling" since 
all scalings give similar results. We see that the scaling is indeed verified for p, J/vP, and 4> and the data points show 
a tendency to be more dispersed for cf>. Interestingly enough, "Diffusive Scaling" leads to the best QF. 



VI. CONCLUSIONS 



In this paper we analysed the scaling properties of deep inelastic observables such as the total cross section, DVCS, 
vector mes on p roduction, and the diffractive cross section, using all available low x data, i.e. from the HI [2fj], ZEUS 
[2lj | . NMCjUj and E665 [23| experiments. In order to be as model independant as possible, we used the quality factor 
method [12j| . allowing to test the validity of scaling variables without any assumption concerning the shape of the 
scaling function. Four scaling variables (HE]), motivated by various approximations of QCD high energy evolution, 
were considered, including the new "Running Coupling II" one derived in [6j . Different versions of those scalings were 
tested, setting the rapidity shift Yq to and the reference scale A to either 0.2 GeV or to 1 GeV, or taking A and Yq 
as free parameters to be fitted. 

We showed that the "Fixed Coupling", "Running Coupling I" and "Running Coupling II" scaling variables lead 
to a good scaling behavior of the total cross section from the F2 data, with similar quality factors. The "Diffusive 
scaling" is disfavoured compared to them. More precisely, the "Running Coupling Ilbis" , with the free parameters Fo 
and A, gives the best scaling behavior in the range 3 GeV 2 < Q 2 < 150 GeV 2 and x < 0.01, with reasonnable values 
of Yq and A extracted from the fit of the quality factor. The "Fixed Coupling" scaling is less deteriorated than the 
"Running Coupling II" scaling when one adds the points with 1 GeV 2 < Q 2 < 3 GeV 2 to the fit, which corresponds 
to the boundary of the validity range of the perturbative scaling predictions. Dividing the data in four Q 2 bins, we 
also studied the Q 2 dependence of the optimal value for the A parameter. "Running Coupling I" gives the most stable 
values and "Diffusive scaling" the less stable, and "Running Coupling II" gives stable values except at low Q 2 . 

The quality factors of the scaling variables were fitted indepcndantly on the DVCS data, and give similar results 
as on the total cross section. In particular, for each scaling variable, the values of the parameter A obtained from the 
fit in the two cases are very close. 

We used the values of the parameters obtained from a fit to F2 data to test the various scaling variables on the 
diffractive cross section and vector meson production data. For the former observable we tested both the fixed f3 
scaling behavior in xp and the fixed xp scaling behavior in (3. At fixed ft, we found a scaling behavior up to ft — 0.65. 
For most ft bins the best results are given by "Running Coupling " scaling variables, and the "Diffusive scaling" is 
often disfavoured. 

At fixed xp, the scaling behavior of the diffractive cross section as a function of ft and Q 2 is far less obvious. This 
is not a surprise, as not enough data is available in the genuine small ft region. A tendency of scaling is however 



observed for the xp — 0.03 bin. 

Concerning p, J , and (f> production, we found a reasonable scaling behavior for all tested scaling variables, with 
the hard scale Q 2 + M v , borrowed from vector mesons wave function studies. Surprisingly, the best scaling is for all 
three vector mesons the "Diffusive scaling" . 

As a phenomenological outlook, it seems useful to work out models dipole amplitude wich could incorporate the 
successful scaling laws following the example of e.g. the Iancu Itakura Munier model [27l |. The quality factor method 
is a good tool to detect the scaling properties of data, and can guide the formulation of models with good x 2 . Our 
study is a good incentive for the formulation of appropriate models with running coupling scaling. There exists 
phenomenological models based on fixed coupling diffusive scaling 28] . The test of scaling and the formulation of 
model will be important to analyse future LHC data. In particular, the wider kinematical range open by the LHC 
could help disentangling the scaling laws in competition, or reveal a new scaling (such as the diffusive one) at higher 
energies. 

On the theoretical ground, our phenomenological analysis can help to improve the theoretical analysis of scaling. 
Indeed, a mere comparison between the values of A and the corresponding prediction based on the leading order 
BFKL kernel shows that one should go beyond this approximation to get an agreement. Thus, our determination of 
the basic scaling parameters can help the theoretical QCD analysis. 
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FIG. 5: F2 data: Normalised QF to 1. as a function of A and scaling curve with A fixed to the value corresponding to the best 
QF for "Running Coupling Ilbis" . A Q 2 > 3 cut was applied to the data. 
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FIG. 6: F2 data: Normalised QF to 1. as a function of A and scaling curve with A fixed to the value corresponding to the best 
QF for "Diffusive Scaling" . A Q 2 > 3 cut was applied to the data. 
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FIG. 7: F2 data: A dependence of the normalised QF to 1. for "Fixed Coupling", "Running Coupling I", "Running Coupling 
II", "Running Coupling Ilbis" and "Diffusive Scaling". A Q 2 > 1 cut was applied to the data. 
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FIG. 8: F2 data: scaling plot for "Fixed Coupling", "Running Coupling I", "Running Coupling II", "Running Coupling Ilbis" 
and "Diffusive Scaling". A Q 2 > I cut was applied to the data to obtain the values of the A parameter. The data with Q 2 < 1, 
not included in the fit, are also shown as grey points. 
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FIG. 9: F-2 data: Comparison of the QF for the different scaling as a function of A. A Q 2 > 3 (left plot) or Q 2 > 1 (right plot) 
cut was applied to the data. 
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FIG. 10: F2 data: Variation of lambda as a function of Q 2 for the different 



scalings. 
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FIG. 11: DVCS data: Scaling curves for DVCS data and for "Fixed Coupling" and "Running Coupling Ilbis" 
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FIG. 12: DVCS data: A dependence of the normalised QF to 1. for DVCS data and for "Fixed Coupling", "Running Coupling 
I", "Running Coupling II" and "Diffusive Scaling". The star indicated the values of A obtained with a fit to F2, Q 2 > 3 GeV 2 . 
The values of A obtained in a fit to DVCS data are similar to te ones obtained in a fit to F2 data. 
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FIG. 13: F 2 " data: Scaling curves obtained for fixed values of (5 for F 2 and fixed coupling. The parameters are fixed to the 
values obtained with a fit to F2 data with Q 2 > 3 GeV 2 . 
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FIG. 14: F 2 D data: Scaling curves obtained for fixed values of xp for F 2 D and "Fixed Coupling" . The parameters are fixed to 
the values obtained with a fit to F2 data with Q 2 > 3 GeV 2 . 
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FIG. 15: Vector meson data: Scaling curves obtained for fixed values of xp and for "Fixed Coupling", for vector meson 
data. The parameters are fixed to the values obtained with a fit to F2 data with Q 2 > 3 GeV 2 . 
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TABLE II: F2 data: QF (multiplied by 10 3 ) and pamameters of the fixed coupling, running coupling I and II, and diffusive 
scalings for the different data sets. We distinguish the data sets for Q 2 > 3 and Q 2 > 1 GeV 2 , and we compare the fit results 
using the full data set, or the HI or ZEUS data only. 
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TABLE III: F2 data: A dependence as a function of Q 2 - The QF are multiplied by 10 3 for simplification. 
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TABLE IV: DVCS data: Values of QF (multiplied by 10 3 ) and fit parameters for the different scalings. The paramters 
obtained for all scaling (except "Running Coupling Ilbis") are close to those found for F2 (see Table II). We note that the small 
amount of data lead to a bad precision of the fit for "Running Coupling Ilbis" when three parameters are used. 
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TABLE V: F 2 data: Values of QF (multiplied by 10 3 ) for the different scalings. The parameters A, Yo and A are fixed to the 
values obtained in the fits to F2 data. Cuts on data xp < 0.01, 5 < Q 2 < 90 were applied. 
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TABLE VI: Vector mesons data: Values of QF (multiplied by 10 3 ) for the different scalings. The parameters A, Yo and A 
are fixed to the values obtained in the fits to F2 data. 
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